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Abstract 

Different treatments for the inclusion of higlier-order perturbative QCD cor- 
rections in parton based transport models of relativistic heavy-ion collisions 
are studied and their influence on experimental observables is investigated. 
At RHIC, particle multiplicities may vary by more than 30%, depending on 
the correction scheme. A detailed analysis of the squared relative momentum 
transfers to be expected in parton (re)scatterings at RHIC casts doubt on 
the applicability of pQCD based transport approaches at the SPS and rules 
out the application of higher-order correction schemes based on a rescaling of 
in the running coupling constant a<j at all energies. 
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One of the primary motivations for the study of ultra-relativistic heavy ion colhsions 
is the expectation that a new form of matter may be created in such reactions. That 
form of matter, called the quark-gluon plasma (QGP [^), is the analogue within nuclear 
matter of the plasma phase of ordinary atomic matter. It is expected that experiments at 
the relativistic heavy-ion collider (RHIC), which is scheduled to begin operations later this 
summer, will observe and identify this new state of matter and engage in systematic studies 
of the properties of strongly interacting matter at high temperatures. 

Transport theory has been among the most successful approaches to the theoretical inves- 
tigation of relativistic heavy-ion collisions. Microscopic transport models aim at describing 
the full time evolution from the initial state of the heavy-ion reaction, i.e. two colliding 
nuclei in their ground states, up to the freeze-out of all particles after the reaction. If this 
description is to contain a quark-gluon phase, as it will be necessary for the RHIC energy 
domain, both partonic and hadronic degrees of freedom must be included. In the partonic 
domain, the Parton Cascade Model (PCM) and other approaches |^,|^ based on per- 
turbative quantum chromodynamics (pQCD) and the Glauber approximation to multiple 
scattering theory have been widely used. The numerical implementations of these models 
I^J^ make use of the hard parton scattering cross sections and radiative emissions as treated 
in the established Pythia event generator |p. 

It is well known that higher-order corrections to leading-order pQCD cross sections are 
large and process dependent . Often these corrections are applied in the form of a so-called 
K-f actor which renormalizes the leading-order pQCD cross section and accounts for higher- 
order corrections |jlO|]. Although K-factors have been calculated for many pQCD processes. 



in many cases they have been determined phenomenologically by comparisons with data. 
The choice of the K-factor is particularly important in the case of parton cascade models 
||7|JTT[|, because it enters into the transport cross section and thus strongly influences the 



number of parton rescatterings. A change in the value of the K-factor does not just lead to 
a multiplicative renormalization of the cross section for an observable. 

Two well-established methods have been employed to implement the K-factor in pQCD 
calculations: 

1. in the form of a constant multiplicative factor: 
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= K X (1) 

with values for K in the order of 1.5 < K < 3 for hard parton-parton scattering. 
In this implementation the K-factor is independent of the virtuality scale of the 
scattering. 

2. by changing the scale in the running coupling constant as in the form 

7 full I LO 

-^iQ') = ^^Mv X g^)) (2) 

with values for rj in the order of 0.03 < r] < 0.08. 



The second method, originally proposed by Ellis et al. [Q, introduces a strong 
Q^— dependence, due to the functional form of the running coupling constant a^. In leading 
order: 

dcTpQCD o . rn 12 TT 



dt = (33 - 2nf) In(gVA^) 

How do the two methods compare to each other? Figure |1] shows the ratio of the hard 
parton-parton scattering cross section corrected by a scale change in Q^, a{riQ'^) and the 
scattering cross section corrected by a constant K-factor, cr{K), as a function of the squared 
momentum transfer Q^. The solid line depicts the ratio with values for rj and K used 
in the parton cascade model (code version VNI 4.12) |0,|ll|. The dotted line shows the 
same ratio with parameter values as suggested in Pythia (code version 5.7)Q. For high 
values of Q^, on the order of 100 GeV^, the cross section ratio converges to unity and both 
correction methods yield the same result. For small values (Q^ < 10 GeV^), however, the 
rescaling method leads to a rapidly increasing cross section that diverges prematurely. 
Therefore, one should consider the lower boundary for the applicability of the Q^-rescaling 
method to be at approximately ~ 100 GeV^, at which the deviation between the two 
correction methods exceeds 50%. 



^The Q^-rescaling method is the default in VNI, while the default in Pythia is K = 1 and 
higher-order corrections are an optional choice. 
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The choice of the correction scheme is always critical in transport theoretical applications, 
because a wide range of values is probed. Figure |^ shows the distribution dN/d^Q"^) 
of "hard" parton-parton scatterings calculated for central Au+Au collisions at RHIC (open 
circles). The distribution for central Pb+Pb collisions at the Super-Proton-Synchrotron 
(SPS) at CERN is also shown (black squares). Both calculations were made with the parton 
cascade model code VNI 4.12, using the constant K-factor option with K = 2.5. The Q^- 
scale is defined as = s for s-channel processes and as = ^{m^^ + 'm'j_2) for t- and 
M-channel processes. At RHIC, the average squared momentum transfer {Q'^)=8.3 GeV^ 
and the distribution exhibits a power-law behavior da/dQ"^ ~ ^/Q^i which is due to the 
Rutherford-like [da/dQ"^ ~ ^/Q'^) form of the leading-order pQCD cross section and the 
parton distribution functions which add additional inverse powers of Q. 

In the limit of small values the distribution is dominated by the cut-off value Qq = 
2.1 GeV. For a constant K-factor the average value of the strong coupling constant as{Q^) 
probed in these hard parton-parton collisions is (a^) = 0.29. If the rescaling method 
is used instead, the number of hard parton-parton scatterings increases strongly, since the 
hard parton-parton cross section increases by a factor of four at {Q"^) ~ 8 GeV^ (see Fig. p. 
The average value of the strong coupling constant probed with the rescaled values of is 
(as) = 0.96. It is therefore highly questionable whether the rescaling method for the 
higher-order pQCD corrections is suitable at RHIC energies. 

The situation at the SPS is even more dramatic: the average squared momentum transfer 
{Q"^) here is only 1.9 GeV^. Again, the distribution is limited for small by the cut-off 
parameter which is now chosen to be Qq = 1.0 GeV. Obviously, the entire distribution in 
the SPS energy regime is dominated by the choice of the cut-off parameter Qq and therefore 
the applicability of pQCD in general and the PCM ansatz in particular is very questionable. 
The average value of the strong coupling constant as probed here is (as) = 0.39. Employing 
the rescaling method at (Q^) = 1.9 GeV^ leads to a hard parton-parton cross section two 
orders of magnitude larger than for the choice K = 2.5, and (ag) probed with the rescaled 
values of increases to 9.3, clearly well beyond the domain of applicability of pQCD. 

The situation with respect to the distribution of values is even more complex if 
one takes the initial non-equilibrium nature of an ultra-relativistic heavy-ion collision into 
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account: initial, parton-parton scatterings are bound to occur at far larger values of 
than rescatterings in the later reaction stages when a thermalized parton gas has been 
formed. Figure § shows the time-evolution of the average squared momentum transfer (Q^) 
for central Au+Au collisions at RHIC. The analysis is based on a VNI 4.12 calculation, 
again utilizing a constant K-factor of = 2.5. Early rescatterings of incident partons occur 
at rather large values of > 10 GeV^. However, between the time tc.m. ~ 0.5 fm/c and 
^c.m. ~ 1-5 fm/c the average (Q^) decreases rapidly and then saturates at a value of about 
7 — 8 GeV^ until hadronization sets in. The time span until saturation, Atc.m. ~ 1 fni/c, is 
in agreement with previous predictions concerning the thermalization time of a deconfined 



parton gas formed in ultra- relativistic heavy-ion collisions |T3|. Most importantly, however, 
the strong time-dependence of (Q^) shows that there is no single scale which can be 
chosen to unambiguously set the pQCD scale in an ultra-relativistic heavy-ion collision. 

How does the choice of the higher order correction scheme influence experimentally ob- 
servable quantities at RHIC? Figure ^ shows the rapidity distribution of negatively charged 
hadrons, calculated with VNI 4.12 and a constant K-factor K = 1.0 (full circles), K = 2.5 
(open squares) and rescaling (full triangles). Introducing higher-order pQCD corrections 
via a constant K-factor K = 2.5 leads to an increase in the h~ rapidity-density at yc.m. of 
roughly 50%, performing a correction based on a rescaling would increase the multiplic- 
ity by an additional 30%. Notably, the changes only affect the area around mid-rapidity, 
—2 < Hem < 2 since the beam- and target-rapidity regions are dominated by soft beam-jet 
physics and are not affected by the amount of hard parton rescattering. Even though the 
difference in the hard parton-parton cross section between the two methods is a factor of 
four at {Q"^) ~ 8 GeV^, the change in final hadron multiplicity amounts to only 30%. This 
"reduction" in the observable effect compared to the change in cross section is caused by the 
interplay of hard scatterings and soft beam/ tar get- jet fragmentation processes. An increase 
in the amount of hard scattering may lead to a decrease in the soft contribution so that the 
sum of both is lower than one may expect from the scaling of the hard cross section alone. 

Based on the above findings, the use of the rescaling method for implementing higher- 
order pQCD corrections should be regarded as unacceptable in the framework of the parton 
cascade model. A major revision of the PCM code VNI, which has employed this method as 



5 



the default, is in progress [0. A forthcoming pubhcation will present a thorough analysis of 
elementary hadron-hadron data for the determination of the appropriate phenomenological 
K-factor, as well as a full set of predictions for nucleus-nucleus collisions at RHIC. One of the 
many intriguing issues to be addressed is the question whether the hot and dense partonic 
matter state created in these collisions introduces a medium-dependence into the K-factor, 
or - more generally speaking - whether it will be possible to probe the medium properties 
of hot and dense QCD-matter at RHIC by means of pQCD processes. 

In summary, we have investigated different higher-order pQCD correction schemes cur- 
rently used in transport models in the context of ultra-relativistic heavy-ion collisions at 
RHIC and SPS. Measurable quantities like the multiplicity of negative hadrons at RHIC 
are influenced on the order of 30%, depending on the choice of the correction scheme. An 
analysis of the regime to be expected at RHIC shows that correction schemes based on 
a rescaling of in the running coupling constant yield unacceptable results due to the 
divergence of this correction scheme for small values of Q^. The application of the parton 
cascade model to nuclear collisions at the CERN SPS which were based on numerical results 



obtained with the VNI implementation of the PCM [jT5[ therefore need to be reconsidered. 

We acknowledge helpful discussions with D.K. Srivastava, K. Eskola, X.N. Wang and 
A. Dumitru. This work has been supported in part by the Alexander von Humboldt Foun- 
dation, by the U.S. Department of Energy (grant DE-FG02-96ER40945), and by the Intel 
Corporation. 
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FIGURES 




(GeV^) 

FIG. 1. Ratio of the hard parton-parton scattering cross section a{r]Q'^) with higher-order 

corrections expressed as a scale change in as{Q'^) (soUd Unc), and the cross section corrected with 
a constant K- factor, a{K) (dotted hne) versus the squared momentum-transfer Q^. The scale 
renormalization leads to a divergent cross section at small, but still reasonable values of Q^. 
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FIG. 2. dN/d{Q'^) distribution of hard parton-parton scatterings in central Au+Au collisions 
at RHIC (open circles) and in central Pb+Pb collisions at the CERN/SPS (full squares). The 
calculations were performed with a constant K-factor {K = 2.5). The distribution shows an 
asymptotic behavior da/dQ^ ~ l/Q^ for large Q^, the distribution at the SPS is strongly cut-off 
dominated. 
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Au+Au @ RHIC 
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FIG. 3. Time-evolution of the average squared momentum transfer (Q^) in hard parton-parton 
interactions for central Au+Au collisions at RHIC. The strong decrease of {Q^) and subsequent 
saturation are an indication of the onset of local thermalization. The strong time-dependence does 
not permit the choice of an unambiguous fixed scale for the entire reaction. 



10 



Au+Au @ RHIC 
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FIG. 4. Rapidity distribution of negatively charged hadrons (h ) from VNI 4.12, calculated 
with K-factors of K = 1.0 (full circles), K = 2.5 (open squares), and rescaling (full triangles). 
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